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Numerical simulation studies are reported for the convection of a supercritical fluid, 3 He, in 
a Rayleigh-Benard cell. The calculations provide the temporal profile AT(t) of the temperature 
drop across the fluid layer. In a previous article, systematic delays in the onset of the convective 
instability in simulations relative to experiments were reported, as seen from the A T(t) profiles. 
They were attributed to the smallness of the noise which is needed to start the instability. Therefore 
i) homogeneous temperature noise and ii) spatial lateral periodic temperature variations in the top 
plate were programmed into the simulations, and A T(t) compared with that of an experiment with 
the same fluid parameters. An effective speed-up in the instability onset was obtained, with the 
best results obtained through the spatial temperature variations with a period of 2 L, close to the 
wavelength of a pair of convections rolls. For a small amplitude of 0.5 /j,K, this perturbation gave 
a semiquantitative agreement with experimental observations. Results for various noise amplitudes 
are presented and discussed in relation to predictions by El Khouri and Carles. 

PACS numbers: 44.25+f, 47.27.Te, 64.70.Fx 


I. INTRODUCTION 

In recent papers, convection experiments of su- 
percritical 3 He at the critical density in a Rayleigh- 
Benard cell with a constant heat current q were 
reported[l, 2], After q is started, the tempera- 
ture drop A T(t) across this highly compressible fluid 
layer increases from zero, an evolution accelerated 
by the “Piston Effect” [3-5]. Assuming that q is 
larger than a critical heat flux necessary to produce 
fluid instability, A T(t) passes over a maximum at 
the time t = t p , which indicates that the fluid is con- 
verting and that plumes have reached the top plate. 
Then truncated or damped oscillations, the latter 
with a period t osc , are observed under certain con- 
ditions before steady-state conditions for convection 
are reached, as described in refs.[l, 2]. The scenario 
of the damped oscillations, and the role of the “pis- 
ton effect” has been described in detail in refs. [6] 
and [7] and will not be repeated here. The height of 
the layer in the RB cell was L = 0.106 cm and the 
aspect ratio T=57. The 3 He convection experiments 
extended over a range of reduced temperatures be- 


tween 5xl0 -4 < e < 0.2, where e = (T—T c )/T c with 
T c = 3.318 K, the critical temperature. The trun- 
cated - or damped oscillations were observed for e > 
0.009 and over this range the fluid compressibility 
varies by a factor of about 30. 

The scaled representation of the characteristic 
times t osc and t p versus the Rayleigh number, and 
the comparison with the results from simulations has 
been described in ref. [8] . Good agreement for the pe- 
riod t osc was reported. However a systematic discrep- 
ancy for the times t p shows that in the simulations 
the development of convection is retarded compared 
to the experiments. This effect increases with de- 
creasing values of [i?a colT — Ra c ] , where Ra corr is the 
Rayleigh number corrected for the adiabatic tem- 
perature gradient as defined in refs.[l, 2] and Ra c is 
the critical Rayleigh number, 1708. This is shown 
in Fig.l of ref. [8], in particular in Fig. lb) for e = 
0.2 and q = 2.16X10" 7 W/cm 2 ([Ra corr - Ra c \ = 
635), where an experimental run is compared with 
simulations for the same parameters. Here clearly 
the profile A T(t) from the simulations shows the 
smooth rise until the steady-state value, AT = qL/X 
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= 95 /i K has been reached, where A is the thermal 
conductivity. Only at t « 90 s. does convection 
develop, as shown by a sudden decrease of A T(t). 
By contrast, the experimental profile shows a much 
earlier development of convection. Fig.l of ref. [8] is 
representative for the observations at low values of 
[Ra corr — Ra c }. At high values, both experiment and 
simulations show the convection development to take 
place at comparable times, as indicated in Fig. 5b) 
of ref. [8], and specifically in Fig. 2 a) of ref. [7], where 
[Ra corl — Ra c ] =4.1 xlO 5 . It is the purpose of this 
report to investigate the origin of this discrepancy 
by further simulation studies. 

II. CONVECTION ONSET 
CALCULATIONS, SIMULATIONS AND 
COMPARISON WITH EXPERIMENTS 

El Khouri and Carles [9] studied theoretically the 
stability limit of a supercritical fluid in a RB cell, 
when subjected to a heat current q started at the 
time t = 0. Their fluid was also 3 He at the crit- 
ical density, and the same parameters as in ref.[l] 
were used. They calculated the time finstab and also 
the corresponding AT(f instab ) for the onset of fluid 
instability and they determined the modes and the 
wave vectors of the perturbations for different sce- 
narios of q and e. For t > finstab inhomogeneities in 
the RB cell and noise within the fluid will produce 
perturbations which will grow, from which the con- 
vection will develop. An indication of the growth 
of convection is a deviation of the A T(t) profile in 
the experiments or in the simulations from the cal- 
culated curve for the stable fluid (see for instance 
Eq.3.3 of ref[6]). It is readily seen from simulation 
profiles such as Fig. la) and b) in ref. [6] that the de- 
viation becomes significant for f only slightly below 
f p - the maximum of AT(t). In simulations, the 
effective start of convection can also be seen from 
snapshots in 2D of the fluid temperature contour 
lines at various times, as shown in Fig. 5 of ref. [10]. 

P.Carles [11] has argued that the reason for the 
discrepancy for the time t p between experiment and 


simulation is that in the former, the physical system 
has noise and inhomogeneities which cause the per- 
turbations beyond f; nstab to grow into the developed 
convection. By contrast simulations have a much 
smaller noise. Therefore in the simulations the per- 
turbations take a longer time to grow than in the 
physical system, leading to a larger t p than observed. 
Carles’ comment led us to try as a first step imposing 
a thermal random noise on the top plate of the RB 
cell, which was to simulate fluctuations in the upper 
plate temperature control of the laboratory exper- 
iment. The temperature of the plate was assumed 
to be uniform, because of the large thermal diffu- 
sivity Dt « 2 x 10 4 cm 2 /s. of the copper plate in 
the experiments. Accordingly simulations were car- 
ried out with a homogeneous time-dependent tem- 
perature random fluctuation of given rms amplitude 
imposed on the upper plate. This implementation 
consisted in adding or subtracting randomly temper- 
ature spikes T t at the time t with a programmed rms 
amplitude at steps separated by 0.02 s. This inter- 
val is much larger than the estimated relaxation time 
of the top plate over a distance 2 L, approximately 
the wavelength of convection roll pair. Values of the 
variance A = yj < (T t — < T t >) 2 > were chosen be- 
tween 0 and 40 /xK. The range of the A values was 
taken well beyond the estimated fluctuation rms am- 
plitude during the experiments [1] of « 1 fj,K/y/ Hz. 
Three representative curves with 0, 3 and 40 /i K are 
shown in Fig. la) by dashed lines for e = 0.2 for q= 
2.16X10" 7 W/cm 2 , L = 0.106cm and T = 5.1. For 
this value of q, the calculation by El Khouri and 
Carles [12] give t insta b = 6.3 s and A T(t instab ) = 75 
/xK. In the simulation without imposed noise, the 
start of convection has therefore been considerably 
delayed relative to fi ns t a b- The injection of random 
noise has a significant effect in developing convec- 
tion at an earlier time. In Fig. la) the three curves 
are also compared with the experimental one, shown 
by a solid line. Here we have not incorporated into 
the simulations the delay affecting the experimental 
temperature recording, so that they could be inter- 
compared more readily, and also with predictions[12] 
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However this operation will be presented in Fig. 4. 
Further simulations with added random noise were 
carried out for e = 0.2 and 0.05 where the A T(t) 
time profiles are not shown here. 

Fig. 2a) shows a plot of the time of the developed 
convection, represented by t v , versus the random 
rms amplitude A for three series of simulations, all 
taken for a cell with T = 5.1. They are a) and b) 
e = 0.2, q = 2.16 and 3.89x 10” 7 W/cm 2 , and c) e 
= 0.05, q = 60 nW/cm 2 , ([Ua corr - Ra c ) = 635, 
1740 and 4200). The simulation results, shown by 
solid circles, are compared with the experimentally 
observed t p shown by horizontally dot-dashed lines. 
It can be clearly seen that noise imposition, which 
creates a vertical disturbance across the fluid layer, 
reduces the time of convection development. While 
the decrease of t p is strong for small values of A, it 
saturates at a certain level of noise amplitude. The 
gap between simulations and experiment increases 
with a decrease of [Ra corr — Ra c ], namely as the fluid 
stability point is approached. A “critical slowing 
down” is seen in the effectiveness of the perturba- 
tions in triggering the instability. Hence this mode 
of noise introduction fails, because its amplitude is 
limited to the vertical z direction and it evidently 
couples only weakly into the convective motion. 

In parallel with the present experiments, S. 
Amiroudine[13] also carried out a systematic study 
of simulations on supercritical 3 He in a RB cell for 
several values of e and q. The resulting profiles 
A T(t) could be compared with those from experi- 
ments done under nearly the same conditions. In 
his simulations, homogeneous temperature random 
noise was again imposed on the top plate. The shift 
in t p showed less systematic trends than in the re- 
sults described in this report. However for the same 
values of e and q as those reported above, and at zero 
noise, the t v values tended to be somewhat smaller 
than in the results of Fig 2a). 

Here we mention that the onset of convection in 
the simulations is further influenced by the aspect 
ratio r. The simulations described above, but with- 
out noise, were carried out in a cell T = 5.1 hav- 



FIG. 1: a) the temperature profile A T(t) from exper- 
iments (solid line with noise) and from several simula- 
tions (dashed lines) at e = 0.2, q= 2.16x10”' W/cm 2 . 
In the simulations, T=5 and uniform temperature noise 
has been imposed on the top plate with variance 
= 0, 3 and 40, as described in the text, b) Tempera- 
ture profile A T(t) from several simulations at e = 0.2, 
q= 2.16x10”' W/cm 2 , T = 8 and imposed lateral pe- 
riodic, time independent temperature variations on the 
top plate with period 2 L and amplitude = 0, 0.5 

and 1.5. 

ing periodic lateral boundaries. Further simulations 
with zero noise for e = 0.2 with T = 8.0, 10.2, 20.5 
and 41.0 were carried out, and showed a decrease of 
the convection development time from ss 90 s, tend- 
ing to a constant value of « 60 s. above T = 20. 
This shift in the onset of instability is probably due 
to the decreased finite size effect which the rising 
plumes experience with increasing T, in spite of the 
periodic boundary conditions. 

The next step in our attempts, stimulated by 
communications with P. Carles, was introducing 
perturbations into the simulations via some time- 
independent lateral variation proportional to sin 
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FIG. 2: a) The time for effective development of convec- 
tion, characterized by t p , versus A (homogeneous tem- 
perature noise imposed on the top plate). The horizon- 
tal dot-dashed lines indicate the observed t p , corrected 
for instrumental recording delay, b) The time for effec- 
tive development of the convection, labeled by t p versus 
B (lateral time-independent periodic temperature vari- 
ations). The horizontal dashed line indicates f p as ob- 
tained by experiment, corrected as before. 

(27 tx/P) where P is the period. A periodic variation 
in the spacing of the plates proved not to be feasi- 
ble for our simulation code, and instead we opted 
to introduce again a temperature variation in the 
top plate with an amplitude B (in /iK) and period 
P=2L, nearly the same as the wavelength of a pair 
of convection rolls. The temperature of the bottom 
plate was kept homogeneous. This “Gedanken Ex- 
periment” implies that the material of the top plate 
permitted a temperature inhomogeneity, which of 
course is not realized in the experiment. As a control 
experiment, we also made a simulation with P = L. 

Fig lb) shows representative profiles A T(t) for the 
parameters e = 0.2 and q= 2.16xlCK 7 W/cm 2 and 
with B = 0, 0.5 and 1.5 /iK, and for T=8. For B = 
0, the effect of having a larger T, namely 8 instead of 
5.1, can be seen by comparing with the curve A=0 
in Fig. la). As B is increased from zero, there is a 
large decrease in the time for convection develop- 
ment, represented by t p , which is plotted versus B 
in Fig 2b). The horizontal dashed line shows the t p 


from the experiment, and this plot is to be compared 
with Fig. 2a). For an inhomogeneity amplitude of 
only B= 0.5/tK, t p is nearly the same for simula- 
tions and experiment. By contrast, simulations with 
B=2/tK and P=L (not presented here) show no dif- 
ference from those with B=0. Hence the nucleation 
of the convection is accelerated if the period is in ap- 
proximate resonance with the wavelength of a con- 
vection roll pair. The values of steady-state AT and 
tosc are only marginally affected by the noise. 

We note from Fig. lb) that the simulation curve 
calculated for B = 0 shows the fluid not convert- 
ing until « 70 s. For the curves with B= 0.5 /iK., 
the start of deviations from the stable fluid curve 
cannot be estimated well from Fig. lb) but is readily 
obtained from the data files, which tabulate A T(t) 
to within 1 nK. For B = 0.5 /iK, systematic devia- 
tions SAT(t,B) = [AT(f,B=0) - AT(t,B)] increase 
rapidly from 1 nK for t >8 s (where AT Ri85 /iK), 
a value comparable with the predicted fi ns tab = 6.3 
s., AT(f instab ) = 75/iK[12]. However a comparison 
with predictions becomes more uncertain as B is in- 
creased. This is because the changes in the base 
heat flow by the periodic perturbations are not con- 
sidered in the theory [12]. We also note that the time 
interval St = [i p — t; n stab] between the first sign of 
instability (SAT > 0) and t p is ~ 20 s, and roughly 
independent of B . This represents approximately the 
period taken by the convection to develop and for the 
plumes to reach the top plate boundary. 

In Fig. 3, we present a series of 2D “snapshots” at 
various times for the simulation with B = 0.5/tK, 
showing the temperature contour lines (in color) for 
the RB cell. The “warm” side is shown by red, 
T(t, 2 = 0) and the “cold” side by mauve, T(z = L ) 
= const. At t= 8 s. the fluid instability has just 
started near the top of the layer, while near the bot- 
tom the temperature contour lines are still horizon- 
tal. At t = 27 s., where the peak of A T(t) at z = L 
has been reached, the warm plumes have reached 
the top plate, and the “cold” piston effect is about 
to start, causing the bulk fluid temperature to drop 
and A T(t) to decrease. The transient process con- 
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FIG. 3: Snapshots in 2D for the RB cell with an aspect 
ratio of 8, of simulations with B = 0.5 q,K at various 
times t after starting the heat current q. The temper- 
ature contour lines and their evolution are described in 
the text. At the time t = 9.9s, the fluid instability has 
just started near the top. 

tinues with damped oscillations of A T(t). Steady 
state convection is reached at t= 80s, with a pair 
of convection rolls having a wavelength of « 2 L, as 
expected. 

In Fig. 4 we show the profiles A T(t) from the ex- 
periment and from the simulations with a periodic 
perturbation amplitude B = 0.5 /iK. For an opti- 
mal comparison, the delay affecting the experimen- 
tal temperature recording was incorporated into the 
simulation curve. For this, the delay function with 
the instrumental time constant r = 1.3 s. [1] was 
folded into the simulation curve by a convolution 
method. This operation retards the initial rise of 
the temperature drop by the order of 2-3 seconds, 
and brings both experiment and simulations into fair 
agreement in the regime where the fluid is stable. 
The time t p for the maximum is now closely the 
same for both experiments and simulations. How- 
ever beyond the predicted instability time finstab= 
6.3 s., the experimental curve starts to deviate more 
rapidly with time than do the numerical simulations 
from the calculated curve for the fluid in the sta- 
ble regime. As discussed previously [2], for these pa- 
rameters of e and q the experiment does not show 
damped oscillations, which are observed for higher 
values of q. In the steady-state, the agreement is 
very good. 
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FIG. 4: Comparison of the profile A T(t) from experi- 
ment and from simulations with B = 0.5 ; iK . To make 
the comparison realistic, the simulations have been con- 
voluted with the same “instrumental” delay time r =1.3 
s. which has influenced the shape of the experimental 
curve. 

Our goal has been to show that injecting a small 
temperature perturbation into the top plate, pro- 
duces for the simulations an earlier start in the con- 
vective instability, which becomes consistent with 
experimental observations. For this, we have lim- 
ited ourselves to an example at a low value of 
[i?a corr — Ra c ], where the delay has been particu- 
larly large with respect to the experiment. 

III. SUMMARY AND CONCLUSION 

We have presented a comparison of numerical sim- 
ulations with experimental data investigating the 
transient to steady convection after the start of a 
heat current through a supercritical 3 He layer in a 
RB cell. Here the temperature drop A T(t) across 
the fluid layer versus time t was studied. The aim 
was to understand and to reduce the discrepancy 
between experiment and simulations in the time of 
the convection development, as detected by A T{t). 
Simulations for one set of fluid parameters (where 
the largest discrepancy had been observed) are re- 
ported with imposed temperature variations on the 
top plate. Satisfactory results were obtained for spa- 
tial lateral temperature variations with an amplitude 
of 0.5 /rK and a period approximately equal to that 
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of the wavelength of a convection roll pair. As the 
perturbation amplitude is further increased, the de- 
velopment of convection occurs earlier than the ob- 
served one. 
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